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The Josephson plasma resonance in layered superconductors with alternating critical current 
densities is investigated in a low perpendicular magnetic field. In the vortex solid phase the current 
densities and the squared bare plasma frequencies decrease linearly with the magnetic field. Taking 
into account the coupling due to charge fluctuations on the layers, we extract from recent optical 
data for SmLai_ a; Sr a; Cu04_ ( 5 the Josephson penetration length A a f, « 1100 A parallel to the layers 
at T = 10 K. 

PACS numbers: 74.25. Gz, 42.25.Gy, 74.72.-h, 74.80. Dm 



The suppression of the Josephson interlayer coupling 
and consequently the Josephson plasma resonance ( JPR) 
frequency, ujq, due to misaligned pancake vortices in lay- 
ered superconductors with identical intrinsic junctions 
was calculated in the Refs. [[!]-§} in good agreement with 
experimental data in the P^S^CaC^Os-a supercon- 
ductor. It was found that in the vortex crystal phase, 
at low fields, uq drops linearly with the vortex concen- 
tration (with the magnetic field B applied along the c- 
axis), while in the vortex liquid and glass state, at high 
fields, ujq drops as 1/ \f~B~ because vortices are uncorre- 
cted along the c-axis and many of them contribute to the 
suppression of the Josephson coupling at a given point. 



Recently the JPR was studied by optical means m 
the layered superconductor SmLai_ x Sr x Cu04_,5 with al- 
ternating intrinsic Josephson junctions containing SmO 
and LaO as insulating barriers between superconduct- 
ing CuC>2-layers. Reflectivity and transmission measure- 
ments were performed with the incidence of light paral- 
lel to the layers (a6-plane) and the polarization of the 
electric field in c-direction. Then the effective dielectric 
function e c ff (co) was extracted from the data by using the 
Fresnel formulae. In crystals with x = 0.2 two peaks in 
the loss function, L(lu) = Im[— l/e c g(uj)], were found at 
frequencies 6.6 and 8.9 cm -1 at B — 0. These frequen- 
cies drop with B linearly at low B < 0.5 T and as \j \f~B 
at B > 1 T. The behavior of the JPR frequencies at high 
magnetic fields was attributed to the vortex liquid state. 

The relative intensities of these peaks and the depen- 
dence of the JPR frequencies on B at high fields were 
explained in the Refs. |lQ-[l2] in the model of two alter- 
nating intrinsic junctions with different Josephson criti- 
cal current densities. Thereby it is essential to take into 
account the charge coupling of neighboring junctions, i.e. 
the c-axis spatial dispersion of Josephson plasmons [j^ , 
in order to explain the peak amplitudes both for paral- 
lel and grazing incidence on the surface [p^[-p3[. The 
dimensionless parameter a characterizing the change of 
the chemical potential of superconducting layers with the 



electron concentration was estimated from the field de- 
pendence of the JPR resonances in the loss functions jll] . 

In this paper we will calculate the behavior of the 
JPR frequencies at low magnetic fields. This al- 
lows us to extract another important parameter of the 
SmLai-ajSr^CuC^-i superconductor, the London pene- 
tration length X a b- 

We consider a crystal with alternating Josephson crit- 
ical current densities J m , m — 1,2, and corresponding 
bare plasma frequencies <J^ m = 87r 2 csJ m /eo ( I > 0: where s 
is the interlayer spacing which we assume to be similar in 
both junctions (for SmLa 1 _. I .Sr.,.Cu04_ ( 5 the difference is 
about 1.5%) and eo is the high frequency dielectric con- 
stant (sa 19 for SmLai-zSr^CuC^-a, see Ref. |11 ). We 
calculate first the dependence of the bare JPR frequen- 
cies LOQ m (B) on B without charge coupling between the 
junctions. Then we will account for the JPR dispersion 
due to a to find the real JPR eigenfrequencies w m and 
the field dependence of the resonances uj m in the loss 
function L(lu). 

Let us consider the field dependence of the JPR fre- 
quencies in the single vortex regime at low fields B <C 
B Jm ,B x , where Bj = $ /A 2 m , B x = $ /47rA 2 b and 
Ajm = ( < E > o cs /87r 2 A 2 b Jm) 1 / 2 are the Josephson lengths 
of the junctions of type m = 1, 2. In this limit the in- 
teraction between the vortices is screened (single vortex 
regime). We neglect pinning assuming that mainly ther- 
mal fluctuations are responsible for the meandering of 
the vortex line. This assumption is valid at sufficiently 
high temperatures T. In Ref. || it was shown that in the 
single vortex regime the field dependence is given as 

u 2 0m (B) = ^ 2 m (0)(l - B/B 0m ), B 0m = $ // m , (1) 
I m = J dr([l - cosy>„,i ;ri+ i_ m ,2(r)]), (2) 

where r is the in-plane coordinate, f n ,i;n-m+i,2( r ) is the 
phase difference between the layers nl and n — m + 1 , 2 
and (...) is the average over thermal disorder. We in- 
troduce the displacements u nm of the vortex from the 
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straight line along the c-axis in the layers nm and its 
Fourier transforms u m (q) with respect to the coordinate 
n of the unit cell consisting of two different junctions, 
< q < 2ir. Then we express the phase difference via 
vortex displacements in quadratic approximation as it 
was done in Ref. 101 and obtain 
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In the following we will use the self-consistent harmonic 
approximation (SCHA) replacing [u n \ — u n 2) 2 by its av- 
erage value ((u n i - u n 2) 2 ) = 4r 2 l5 where r wrn is the 
meandering length for junctions of the type m. I m in 
terms of r wm is 



I m = (7r/2)r 2 m ln(0.8Aj m /r. li , m ). 
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The increase of the vortex energy due to displacements 
is 
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[Ej l \u l {q)~u 2 {q)Y 
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E J2 \u 1 (q) - u 2 (qyo\ 2 + WMQu^q)] 2 + \u 2 (q)\ 2 )}, 

where Ej rn — <f>oJ m /2c is the Josephson coupling den- 
sity in a junction of the type m multiplied by the factor 
(-71-/2) and Wm = ^o 5 / (^^Ib) 2 accounts for the cage po- 
tential due to the nonlocal magnetic interaction between 
pancakes in different layers. After diagonalization of this 
energy we obtain the free energy functional of the vortex 
line 

?= \ E mq)\vM)\ 2 (7) 

g,i=l,2 

E xa {q) = E JX + Ej 2 + W M ± (8) 
(E 2 n + E 2 n + 2E n E J2 cosq) 1 ' 2 . 

After minimization with respect to the new variables 
Vi(q) we obtain the free energy as 

F(E Jm ,W M ,T) = 2Tj2HEi(q)E 2 (q)]+F (T), (9) 
i 

For I m we derive (neglecting logarithmic factors) 
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and similar for I 2 , where £j m — Ej m /WM = 
(A a6 /A Jm ) 2 = (A 2 b /A cm s) 2 and £ = {£ n + £ J2 )/2. The 
values A cm are determined by the zero field bare plasma 
frequencies ujo m (B = 0) = c/y^Xcm in the absence of 
an external magnetic field B. 

For high fields B 3> Bj m , B\ the vortex lattice is more 
dense, the interaction between the vortices diminishes 
pancake fluctuations, while the enhanced tilt stiffness 
of the lattice favours larger fluctuations, i.e. a reduced 
Josephson energy. We write down the energy functional 
in the Fourier representation k with respect to the in- 
plane coordinate r and the unit cell index n. We sepa- 
rate transverse, Ut(k, n), and longitudinal, u/(k, n), vor- 
tex lattice displacements (for details see Ref. @J) 

£vor = £t{u t ) + £l{lli), (11) 

£t(u t ) = [C 66 £: 2 + $44(k,m)]Mk,2n + m)| 2 , (12) 

k , n , m 

£l(ui)= E t* 11 +*4 4 (k,m)]|w ; (k,2n + m)| 2 , (13) 



where the summation over k is limited to the first 
Brillouin zone, which we approximate by the cir- 
cle k < K Q , Kl = 4nB/$ - Further, $ u w 
(B 2 s/47rA 2 b )(l — fc 2 /4if 2 ) is the compression stiffness, 
Ce6 — AQQB$os/(8TrA a b) 2 is the shear modulus. The 
parameter Aqq < 1 describes the fluctuation suppres- 
sion of 6*66, cf. Ref. [l4| , which we approximate as 
= 1 — OAB / B nlc \ t with i? me it being the melting field, 

and 



$44 (k, m) = E cm + E Jm rj m (14) 
is the tilt stiffness. The cage energy is given as 



B<P Q s ( . 0.13a 2 
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Here a 2 = $o/B is the intervortex distance and 



B , 0.11a 2 
r\ m = In ■ 



4tt 



2$ r 2 wm (l-0.bZk 2 /K 2 ) 2 a 4 fc 2 ' 



(15) 



(16) 



After diagonalization of the energy functional we ob- 
tain the eigenvalues E t \ and E t2 for the transverse dis- 
placements with 

E t \E t2 = S t iS t2 + (Sn + St^iEjxrji + E J2 r/ 2 ) + 

2Ej 1 Ej 2 r ll r} 2 (l-cosq), (17) 
Stm = C§§k 2 + E cm . (18) 

For longitudinal eigenvalues Ei m we need to replace Stm 
with S im = $11 + E cm . 

Finally we find the free energy and by differentiating 
it with respect to Ej m we find the self-consistency equa- 
tions for r wm as 
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I m = nr wm /2 = Tr(r wtm + r wlm )/2, 
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The Eqs. (N)-(M) should 



D t = S t \St2 + (S t i + S t 2)(E J1 r) 1 + Ej2m) 
2E J1 E J2 T] 1 ri 2 , 

and similar for r~ t 2 and r 2 ui 
be solved self-consistently to find the meandering lengths 

The bare frequencies u)Q m (B) determine the resonance 
frequencies u> m (B), which are renormalized due to the 
charge coupling of the layers, and are observed ex- 
perimentally in reflectivity transmissivity and in the 
loss function L(u>). Next we derive the renormalized 
frequencies and compare them with the measurement 
of these quantities in Ref. ||. The parameter a = 
(eo/47res)(9/i/(9p) characterizes the interlayer coupling 
due to charge fluctuations where \i and p are the 

chemical potential and charge density on the layers, re- 
spectively. The parameter a was estimated as 0.4 for 
SmLai-^Sr^CuCU-s at x = 0.2 in Ref. both from 
the magnetic field dependence of the JPR resonances in 
the vortex liquid state and from the relative amplitude 
of the resonances in L(ui) for B = 0. For incident light 
perpendicular to the crystal surface ac the wave vector, 
k x , of the wave propagating into the crystal determines 
the reflection coefficient, R(oj), according to the Fresnel 
formula, R = | (1 — n)/(l + n)\ 2 , with the refraction index 
n(u>) = ck x (uj)/uj. In Ref. [^ljjlj] the dispersion relation 
k x (uj) was calculated using the Maxwell equations and 
the equation for the phase differences in a stack of in- 
trinsic Josephson junctions coupled inductively and due 
to charge variations inside the layers. The result is 
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are the normalized JPR frequencies, 
w = u)' 2 /u>qi, cr m = Awa m /eou>oi, and a m are the c-axis 
quasiparticle conductivities of the junctions. We see that 
the loss function L{u>) has two peaks at uj = oj m corre- 
sponding to zeros of e e g (uj) in the absence of dissipation. 
These resonances in the loss function correspond to the 
transverse plasma modes propagating along the layers. 
There is another characteristic frequency defined by the 
relation 



co 2 polc (B) = [uj 2 01 (B)+uj 2 2 (B)}(2a + l/2), (24) 

which corresponds to a pole of Im(e e gf(w)) in the absence 
of dissipation, cf. Eq.(po|). This frequency is near the 
minimum of the loss function L(u>). This frequency is 
often called "transverse plasmon" and it is well defined 
experimentally as the peak in the real part of the optical 
conductivity, Im(o;eeff), see Ref. ||. 

Next we calculate the field dependence of the frequency 
Wpoic, expressed by Eq. (24) via the field dependence 
of the bare frequencies uio mi Eq. ([!]), which we already 
found. We obtain 

^ ole (S)=^ ole (0)(l-5/B p ). (25) 

In the single vortex regime (B <C Bj m , B\) the result is 

= ggeg(cjgi(0) +^o 2 2 (0)) 
p 32n3c 2 sTg(£j 1 ,£ J2 ) ' 
g(£ji,£j 2 ) = (l/2)[f(Ej 1 ) + f(Ej2)] 
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FIG. 1. Field scale B P (B) in the low (B < B Jm ,B x ) or 
high (B » Bj m ,B\) field limit using Eq. EQ) (dashed) or 
Eq. M) (solid, B me i t = 2 T) for \ ab = 1700 T = 2.3 K. 
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FIG. 2. Field scales Bo m and B p for the variation of the 
JPR resonances u m and the pole o; po i c in the limit of high 
fields B > B x , cf. Eq. @ (B mclt = 2 T, A a6 = 1700 A, 
T = 2.3 K). 

In Fig. |l| it is seen that in the low field limit B <C Bj m 
as determined in Eq. ( |i"(i| ) the field scale -B p determining 
the variation of the peak frequency uj po \ c in Im(cje c ff) is 
field independent. In contrast to this, the self-consistent 
solution of Eq. (|l^) valid for B 3> B x introduces a vari- 
ation of B p with B. In Fig. || the field scales B 0m and 
Bp responsible for the variation of the JPR peaks uj m in 
L(u>) and of the pole a; po i in Im(we e ff) are compared. 
The sensitivity of the magnetic field scale Bo m to the 
choice of the correct melting field is low. 

1 t. . . , . . , 



We now compare the calculated frequencies u> m (B) 
with the experimental results in Ref. Q , where the eigen- 
frequencies LO 0rn (B) are extracted from the zeros of the 
effective dielectric function e e g and the characteristic fre- 
quency Wpoio from the peak in Im(weeff). 

Firstly, we use the normalized values of 
Wpoie(-B)/w P oie(0), as they do not depend on the choice of 
the charge coupling parameter a, see Figs. || and ^. Fol- 
low fields the results derived in the cases B -C B\ , Bj m 
(Eq. ©) and B > B x , B Jm (Eq. ©) are close, but de- 
viate at fields approaching the melting field -B me it ~ 0.5 
T, which is chosen in accordance with the observed phase 
transition to the vortex liquid. The agreement with the 
experimental data is worsening near the vortex solid- 
liquid transition. 
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For low temperatures, cf. Fig. ||, agreement with the 
experimental data can be obtained for A a & = 1850 A, 
while at high T = 10 K we get A afc = 1600 A us- 
ing Eq. djj or A afc = 1100 A (B m = 0.5 T) from 
Eq. (|l9|). Following from estimates for B\ the latter case 
seems to be more adequate. For higher temperatures, cf. 
Fig. ||, the assumptions of our theory are expected to be 
better justified and in this case the obtained value for 
X a b = 1100 A is comparable to the ones reported for 
other layered superconductors. 



FIG. 3. Frequency £J po i c /cj po i c (0) of the peak in Im(a;e e ff) 
in the low and high field limit B <C Bj m , B\ or5 > Bj m , B\ 
using Eq. ^ or Eq. @ respectively (T = 2.3 K). Best fit of 
the experimental data in Ref. [6] is generally obtained for the 
penetration depth A a 6 = 1850 A, except for the case, when a 
too large melting field B mc i t = 2 T is chosen for comparison, 
which requires X ab = 2230 A. 
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FIG. 4. Best fit of w P oie/w P oie(0) for high temperature 
T = 10 K using the low and high field limits in Eq. ( |1C| ) 
{\ ab = 1600 A) and Eq. @ (X ab = 1100 A, B mclt = 0.5 Tj! 
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FIG. 5. JPR peaks u m (B) in L(lu 
Im(oj£ G ff) using the low field limit (\ a b ■ 
K, a = 0.4, u cl (B = 0)/c = 6.6 cm 
cm -1 , cf. Ref. [6,15]). 
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FIG. 6. Best fit of w 0m (B)/w 0m (0) from Ref. [6] for high 
temperatures T = 10K using the low (dashed, X a b = 1600A) 
and high (solid, \ a b = llOOA, _B mc it = 0.5T) field approxima- 
tions, Eqs. ( |l(i| ) and ( M ) respectively (a = 0.4). 

Assuming the charge coupling a « 0.4 as estimated in 
Ref. p2| from the vortex liquid state and the shape of 
the loss function for B = 0, the resonance frequencies 
oj m (B,a) in L(w) and cj po i c are shown in Fig. ||. The 
experimental data for the normalized JPR resonances 
w m (5)/w m (0) for high temperatures can be fitted with 
the same choice of parameters as for uj po i c , cf. Fig. |^. 

To conclude, thermal fluctuations in the vortex solid 
state describe satisfactorily the dependence of the 
plasma frequencies in the Josephson coupled system 
SmLai_ a; Sr a ;Cu04_5 with two different alternating junc- 
tions on a perpendicular magnetic field. From the com- 
parison of our theoretical results and the experimental 
data we estimate the in-plane London penetration length 
X a b ~ 1100 A at high temperatures T = 10 K, which is 
similar to other layered superconductors. 
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